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Abstract — We analyze the performance of coherent impulsive- 
radio (IR) ultra-wideband (UWB) channel in presence of the in- 
terference generated by concurrent transmissions of the systems 
with the same impulsive radio. We derive a novel algorithm, 
using Monte-Carlo method, to calculate a lower bound on the 
rate that can be achieved using maximum-likelihood estimator. 
Using this bound we show that such a channel is very robust to 
interference, in contrast to the nearest-neighbor detector. 

I. Introduction 

In this work we consider an impulse-radio ultra-wideband 
(IR-UWB) system. This type of communication is being 
actively developed because of its promising features. It gen- 
erates low-power transmissions over a large bandwidth, hence 
achieves high bit rates with low power consumption. Other 
advantages of IR-UWB are low cost, multi-path immunity and 
precise ranging capabilities. The IR-UWB physical layer is 
adopted as a standard for 802.15.4a personal-area networks 
(PAN) in the 3-5GHz band. 

The current IR-UWB receivers are designed to work with 
Gaussian white noise, and they use nearest-neighbor decod- 
ing which, in case of Gaussian white noise, corresponds to 
maximum-likelihood decoding (MLE). However, in a network- 
ing environment, where there are several concurrent transmis- 
sions of IR-UWB devices, the interference is very impulsive 
and nearest-neighbor decoding is suboptimal. Interference may 
occur because of several competing piconets, or during a 
random access phase in the same piconet. 

We are interested in calculating the achievable rates of a 
coherent IR-UWB channel in presence of white noise and 
impulsive interference generated by a network of IR-UWB 
transmitters. Instead of nearest-neighbor, we focus on the 
optimal, MLE decoder. Although we cannot calculate the 
maximum achievable rate explicitly, we give a novel lower 
bounds on the achievable rate. We also present a simple upper 
bound and we use these bounds to evaluate the performance 
of the channel. 

We further consider an example of a channel with a single 
impulsive interferer. We show that if one uses an MLE 
decoder, in many cases it can mitigate the effect of impulsive 
interference, especially when the interference is strong. This 
is in contrast to IR-UWB receivers with nearest-neighbor 
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decoding which are highly penalized by strong impulsive 
interference (TJ. 

In Section [TT] we present related work. In Section [III] define 
system assumptions. We derive the upper and lower bounds in 
Section IIV1 we illustrate their performance in Section [V] and 
we conclude. 

II. Related Work 

The performance of a channel with additive, non-Gaussian 
noise and nearest neighbor decoder is discussed in fl2|. In 
particular, the case with impulsive interference is discussed 
in 0, [1| where the authors show that the Gaussian approxi- 
mation of the interference is not correct and give a numerical 
model to evaluate the bit-error rate in homogeneous settings. A 
similar problem is found in [4], where fast and more efficient 
methods are developed to handle heterogeneous cases and with 
multipath channels. However, in all of these works, the author 
consider the nearest neighbor decoder, which is not the optimal 
one when the interference is not Gaussian. 

Our work is inspired by [5|, which calculates an upper 
and a lower bound on achievable rates of a non-coherent 
IR-UWB channel. This work is extended in ||6l to evaluate 
the performance of different non-optimal detectors. A similar 
analysis is done in Q for transmitted-reference IR-UWB 
radio. Nevertheless, in all of these works the authors do not 
consider the effects of impulsive interference from concurrent 
transmission of the same type of radios. An appoximate imple- 
mentation of MLE decoder for multi-user interference channel 
is given in (8). However, to the best of our knowledge, we are 
the first to numerically calculate bounds on the achievable rate 
of such a channel. 

III. System Assumptions 

A. Channel model 

We consider a set of 7 + 1 nodes. Node 1 communicates 
with node while the other I — 1 nodes generate impulsive 
interference from concurrent transmissions. All nodes use 
IR-UWB for communication, meaning that they transmit a 
sequence of very short pulses. All signals and channels are 
assumed to be real-valued. The received signal at node is 
equal to 

oo 

r(t) = ^2 ui[n]A x hi{t - nT s ) + z(t) 

n— — oo 
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+ Ui[n]Aihi(t-nT s + A in ), (1) 
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where Uj„ is the n-th transmitted symbol by node i, Aj 
is the received amplitude of a pulse at node 0, transmitted 
from node i, A,„ is the time shift (due to asynchronicity) 
between interferer i and destination for the n-th symbol, 
T s is the symbol duration, and hi is the normalized channel 
impulse response of a channel from node i to node 0, which 
is described next. 

When a transmitted signal propagates from one point to 
another, it travels over multiple paths. If a perfect impulse 
5(t) is transmitted, the received signal will be 

L 

h*(t) = ^a,8(t-T a ), 

s=l 

where a s is the attenuation and t s is the delay of s-th path (s 
taking values from 1 to L). We assume that hi(t) is normalized 
so that^fL, a 2 s = 1. 

The received impulse response h p (t) is further filtered to 
the system bandwidth W, and the received signal is 

L 

h(t) = V2W^2a s smc2W(t - r s ). 

s=l 

We further assume that we sample the channel impulse 
response at 2W, and we have the following samples 

L 

hm = h \2w) = V2^E a - sinc2 ^(^7- T -)' 

s—l 

for m = 1, • • • , M where M is the number of samples. Since 
the number of paths is very large, h m is a sum of a large 
number of random variables thus it can be assumed Gaussian. 
However, different samples are not independent. Therefore 
we assume that each h; = (ha, ■ ■ ■ ,hiM) is a multi-variate 
Gaussians and {h^}; are i.i.d. with zero mean and covariance 
matrix T^. 

Channels h; are independent but their statistics depend on 
the time jitter Ai„ = 0. In an unlikely case that interferer i is 
symbol-level synchronized to a receiver (that is A J?l = 0), the 
receiver will receive the full energy of this interferer. Then, 
the distribution of will be the same as the distribution of 
hi, for the intended transmitter, that is, perfectly synchronized 
as well. Symbol-level synchronized interferer is thus a worst 
case approximation. Although we can readily use our model 
for arbitrary A in , we are interested in deriving a lower-bound, 
hence we will assume A in — for all i,n. Hence, we will 
assume that all channel responses have the same covariance 
matrix T, = T. 

For each transmitted symbol we have M channel samples. 
Each sample can be interpreted as one dimension of a channel. 
We can thus formulate the corresponding channel as discrete 
vector channel, and we have the following channel model 

i+i 

r[n] = ui[n)Ai[n]hi[n] + Uj [n] Aj [n] hj [n] +z[n]. (2) 

where hj[n] are the samples of channel impulse responses 
during transmission of symbol n and z[n] are samples of white 
noise hence i.i.d. Gaussian. 



The typical channel coherence time for a UWB channel is of 
order of tenths of milliseconds J9). In practice, this means that 
the channel is constant throughout the duration of a packet and, 
in Equation (O, h^[n] = are constant in n. Furthermore, 
since we consider coherent communication, we assume the 
receiver knows transmitter's channel hi but it does not know 
the interferer s' channels {hijj^i. 

B. Transmitter and Receiver Structure 

It has been shown that an efficient IR-UWB system needs 
to transmit infrequent pulses [ 1 1 due to a very low available 
average transmission power. This in turn facilitates multi-user 
communication. One implementation of such a principle are 
for example time-hopping (TH) codes IfTTI . 

We use a generic signalling model and assume that {ui[n]}„ 
are i.i.d. random variables with P(ui[n] = 1) = and 
P(ui[n] = 0) = 1 — 77,;. Variable rji thus denotes the duty-cycle 
of node i. This model is more general than the TH model as 
it does not impose any dependency between symbols, which 
is needed in case of TH due to implementational constraints. 

We also assumed that the transmitted symbol can be only 
or 1. This is not the most general type of modulation. For 
example, having Ui[n] S {—1, 0, 1} would allow us to transmit 
additional information in the symbol phase. However, for the 
simplicity of the presentation we assume Ui[n] € {0, 1} and 
we note that the result can easily be extended to different 
modulation schemes. 

The average power during transmission is upper-bounded by 
pMAX ^ w i! ere tjjg value of P^ax j s S p ec ifi ec j by regulations. 

Since the probability of transmitting 1 is i^, the amplitude of 
a pulse is bounded by A 2 < Pf 1AX jr\i. 

We shall not assume any particular coding or detection 
techniques for our system. We will use random coding to 
derive a lower bound and Shannon capacity to derive an upper 
bound on the achievable rates. The precise descriptions are 
given next. 

IV. Bounds on Achievable Rate 

We want to derive bound on the achievable link's physical 
data rate, given the received signal power and the received 
powers of interferers. For that matter, we will consider the 
discrete vector channel model described in @. We will first 
present a simple upper-bound, and then we will derive a novel 
lower-bound, which is the main result of our paper. 

A. Upper-bound 

To derive an upper-bound we consider the information- 
theoretic capacity of the channel, constrained to the binary 
input alphabet. We can represent our channel as an M- 
dimensional vector channel 

R=U 1 A 1 h 1 + J2 U i A i h i + z > ( 3 ) 

i>l 

where {Ui]i = i...i are i.i.d Bernouli random variables with 
P(Ui = 1) = rji, z = {z m } m —i...M are i.i.d Gaussian with 
variance crj^, {.1,}, ;.../. {h,„, }, ,.../.„, ,.... u are known 
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signal attenuations and unknown (except for hi) but constant 
channel fadings from transmitter i to the destination respec- 
tively. 

We upper-bound the capacity of the channel assuming 
the receiver knows the received symbols {u;}i>i from the 
interferers. Since it can then perfectly estimate the interference 
and extract it, the only remaining noise is the white noise z. 

We use the notation H = {hi}i = i .... j and U_i = 
{Ui}i—2,— ,i- We can write the upper-bound as 



C u = /(ui;R|hi) 



E 



hi,£/i,z 



M 



log [ P(Ui^v)e 1 

,i>={0.1} 



(4) 
(5) 

(6) 



The result is similar to the result from [6] for a non-coherent 
channel. This bound can be easily calculated using e.g. Monte- 
Carlo simulations. 



B. Lower-bound 

1) Threshold Decoding: Next, we will derive a lower- 
bound on achievable rates using a practical decoding scheme. 
We suppose the source sends data in packets of length TV, 
where N is assume large. Each packet has a coding rate 
Cr associated with it, yielding error probability of decoding 
P(err). We consider an upper bound on P(err) using random 
coding bound technique. 

Suppose packets are coded using a random codebook C C 
{0,1}^, ||C|| = 2 NCr . A source, knowing the channel state 
hi, sends codeword Ui and a destination receives R = 
{r[?i]} Tl= i.... jv, as described in (|2]). The optimal decoder is 
the maximum likelihood (MLE) decoder which selects the 
codeword Ui that maximizes the likelihood P(R|ui,hi). 
However, the performance of the maximum likelihood decoder 
is hard to analyze. Since we are interested in an upper-bound 
on the probability of error, we shall consider a simple threshold 
decoding scheme, based on an arbitrary threshold 0. If the 
likelihood P(R|ui,hi) > for only Ui e C, then the 
decoding is successful. Otherwise, it fails. 

2) Performance Analysis: We start by giving all the no- 
tation we will be using: H = {hj}j=i.... ./, Hi = 
{hi}i =2 , ...,/, R = {r[n]} n =i,...,jv. u, = {ui[n]} n= i ) ... l N, 
U = {Ui}i=i,...,/, U_i = {uf}i =2 ,...,/ and U[n] = 
{ui[n]}i = i,... j. We will also use a short notation for 
P(Y|v,hj = P(R = Y|ui = v,h! = h) and 
P(Y | V, h) = P(R = Y | U = V, hi = h). 

We first need to choose the threshold which will yield 
good performance. Ideally, 0(v, hi) is a function of the 
(unknown) transmitted codeword v and a channel-state hi, 
and we shall choose it to minimize the probability of false- 
negative P(P(R|v,hi) < 0(v,hi)|ui = v,hi). We will 
show later that the optimal 9 does not depend on the choice 
of v, hi. 



The noise and the interferences are ergodic processes hence 
for a large packet size N we have that P(P(R|v, hi) > 
0(v,hi)|ui=v,hi)^lif 

0(v,hi) = (l-e)E R (P(R|v,hi)|ui =v,hi) (8) 

= (l-e)J P(R = Y|v,hi) 2 dY (9) 

for any e > 0. We will choose 0(v, hi) = / Y P(R = 
Y|v,hi) 2 dY (i.e. e = 



0), and assume further 
P(P(R|v,hi) > 0(v,hi)|ui = v,hi) = 1. We next 
show that 0(v, hi) does not actually depend on v, hi, hence 
we can write 0(v, hi) = 0. 

Proposition 1: The following integral 

p(||v-w||, hi) = j P(R = Y | v, hi)P(R = Y | w, hi)dY 

depends only on ||v — w||. Also, 0(v, hi) depends neither on 
v nor on hi. 

Proof: Let us denote with Q[n] = 

Ei=2 u iN^N h iN + Z W- Then ' ^( R - Y|v,hi) = 
P(U n =i-N QN = Y W - v[n]Aihi) and 

J P(R = Y| v,hi)P(R = Y| w,hi)dY = 

p( |J Q[n]=Y[n]J x 

\n=l"JV / 

|J Q[n] = Y[n] - (w[n] - dy. 

The distribution of the vector {Q[n]} n= i...jv is by definition 
symmetric and invariant to a permutation of its elements, 
hence the value of the integral depends only on ||v — w||. 
Furthermore, if | |v — w| | = 0, as in ©, then the integral does 
not depend on hi either. ■ 
Now we are interested in the probability of error of decod- 
ing a random transmitted codeword. We consider a random 
codebook C, and from there select a random codeword v to 
transmit. Note that P(v = o> | C = C) = 2- CrN 1{lj e 
C} since all the codewords from C are equiprobable. The 
probability of error can be bounded by the union bound as 

P(err|hi) < (10) 

< E c ,vec 



P 



E 



P(P(R|«,hi) >0|v,hi) 



(11) 



= 2 c « JV E v , w , v#w [P(P(R|w,h 1 ) > 0|v,hi)] (12) 

where v, w are two randomly choosen codewords from a 
random codebook. 

Next, using Markov inequality, we bound 



P(P(R|w,hi) >0|ui = v,hi)< 
^E R [P(R|w,hi)|v,hi] 
p(||v-w||,hi) 



< 



(13) 
(14) 

(15) 

where the last equation follows from Proposition Q] The 
Markov bound is the best bound we can use knowing only the 
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mean of a random variable, and numerical results in Section [V] 
show that the bound is useful for performance evaluation of 
the channel. 

Random codewords v, w can be assumed independent in 
a large codebook (when CrN is large). We then have that 

P(||v - v 

m) 2 ) N - d - 

Next, let e<i be a vector with d ones and n — d zeros and 
let us denote J(V,W,hi) = / y P(R = Y|V,hi)P(R = 
Y | W, hi)dY. Then from G2) and Q3) we have 

= E[J([e d ,U_i] 



=3.00 



jcuuuk ^wiicii ft±\ is laigc^. vvc men nave liicil 

r\\ = d) = ( ^) (2^(1 - m )f{rjt + (1 - 



p(d) 
P(err) 



2 c R 



N 



N 



(16) 
(17) 



[e , V_i],hi)], 
< ^- r 2^P(\\u-v\\ = d)p(d) 

where U_i = {«i}i=2,-,7, V_i = {«j}i=a,...,/ are ran- 
dom codewords transmitted by interferers. We can express 
J(U,V, hi) in a closed-form, as explained in Appendix, 
and calculate the mean using Monte-Carlo simulations. Since 
p(0, hi) = 8, we can use the same procedure to calculate 8 
(note that in addition 8 does not depend on hi). Details on 
Monte-Carlo simulations are given in Section IV- Al 

From (fT7b we can obtain a lower-bound Cj on the commu- 
nication rate Cr. When N — > oo we can obtain arbitrarily 
small P(err) using communication rate 



C, = --log 2 



£P(||u-v||=d) 



p(d) 



(18) 



\d=0 



V. Numerical Results 

In this section we first discuss the reliability of the Monte- 
Carlo simulation, and then illustrate the results on a channel 
with a single interfere! - . 



A. Monte-Carlo Simulations 

We calculate 8 and ^2^ =0 P(\\u — v|| = d)p(d) using 



Monte-Carlo simulations, averaging over many random sam- 
ples of V_i,U_i and hi. 

In the case of 8, we verify that the samples obtained by 
Monte-Carlo fit the Gaussian distribution well. This allows us 
to calculate confidence intervals of the simulation [12], and 
in all cases the relative confidence intervals are smaller than 
10%. 

In case of ^2^ =0 P(\\u — v j j = d)p(d) the samples are no 
longer Gaussian, but we verify that a log transform is Gaus- 
sian. There is a simple intuitive explanation for this. For very 
small d(d = 1,2) the candidate and the transmitted codewords 
are similar, the probability of error (estimated through p(d)) 
is high. However, there are a few such codewords. On the 
contrary, for large d, p(d) is small, but there are a lot of such 
words. 

In all the simulations we find that the relative confidence 
intervals for the error probability of decoding P(err) are 
smaller than 50%. We are interested in C\ ( TT~8b which is 
of order of log 2 (P(err)). Since the values of interest of 
log 2 (P(err)) are smaller than -20, we can see that the relative 
confidence for C; is approximately 5%. 
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Fig. 1 

Upper and lower bounds on achievable rates of channels with 
a single interferer for i = 3m, i = 5m, and i = 8m. the lower 
bound is given for different values of t]2 (the upper bound is 
independent of r)2>. the upper bound coincides with the 

LOWER-BOUND FOR d = 10UM. 



We also find that the lower-bound coincides with the upper- 
bound in the case of purely Gaussian interference (since our 
lower bound then coincides with a well-known random-coding 
bound for AWGN channels). 

B. Single Interferer Channel 

In order to illustrate our results, we consider a channel with 
a single interferer (/ = 2). We take N — 80 which is around 
the maximum packet size we can use to have reasonable 
simulation times. We also take M = 5 corresponding to a 
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5 -tap receiver, which is also aroung the maximum M we can 
simulate. 

For channel statistics we use the measurements from |6| 
which says that the tap energy drops linearly with the tap 
delay. Approximately 14% of the total energy is in the first 5 
taps. We also assume here that taps are independent. 

In this section we are interested in comparing performance 
of channels for different channel parameters. In order to avoid 
additional unnecessary variance in our results, we will assume 
here that hi is given. Similar results hold for different values 
of hi (hence will hold for the average channel realisations as 
well). 

We further suppose that the distance between receiver 1 
and receiver is I and the distance between interferer 2 and 
receiver 1 is d. Both 1 and 2 transmit with transmitting power 
^trans = 0.1 mW. The average received power at distance I 
is Prcv(0 = ftrans&J-" where b = 10 _5 - 5 ,a = 3.3 are 
taken from 1131 . White-noise power is = lCP 13 !^. The 
maximum communication range of such a system is around 
10m, so in our simulations we use link sizes of I =3,5, 8m. 

In order to illustrate the result we consider a scenario in 
which we fix rji — 0.5 and for different I we vary d and 
r]2- The results are shown in Figure Q] We see that for I = 
3m the achievable rate is insensitive to the interference. For 
I — 5m and / = 8m the rate decreases when the interferer is 
closer than d < 30m. However, it drops until d = 5m, and 
for smaller d the achievable rate stays constant. This shows 
us that when the interferer is close enough, MLE decoding 
actually performs a kind of multi-user detection, successfully 
extracting the interference and preventing a further drop in 
performance. 

We further compare the achievable rate on a channel with 
interferer to the achievable rate without interferer (or equiva- 
lently sufficiently large d, e.g. d — 100m). We see that even 
when the interferer is much closer than the transmitter, the 
rate drop is not significant. In our examples, it is never more 
than 50%. This is in contrast to [fl] where the rate drops to 
zero if interference is much stronger than the signal itself. 

VI. Conclusions and Future Work 

We analyzed the performance of coherent IR-UWB chan- 
nel with MLE detector. We presented a novel procedure to 
calculate a lower-bound on achievable rates using random- 
coding techniques and Monte-Carlo simulations. Using this 
bound we are able to show that the performance of an MLE 
detector is significantly better than the performance of a 
widely used nearest-neighbor detector in presence of a strong 
impulsive interference. This suggested that the use of a more 
complex MLE detector may completely eliminate the need for 
a medium access protocol in IR-UWB networks. The analysis 
of more complex networking scenarios and different medium 
access protocols is left for future work. 
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Appendix 

In the appendix we explain how to calculate the channel 
output distribution P(R = Y|V,hi) and J(V,W,hi) = 
JyP(R = Y|V,hi)P(R = Y|W,hi)dY. First, condi- 
tional to the channel realisation H and the transmitted symbols 
U, the channel outputs R = {r m [w]} are Gaussian i.i.d. RV 
with distribution 

/ v MN 

P(R = Y|U,H) ' 1 



w , 



(N M I \ 

w — 1 m— 1 i— 1 / 

Also, each channel response h^ is multivariate Gaussian with 
distribution 

/ x 7-1 

P(H_i) = ' 



cxp (-i^hfT-O, 



Thus, we have 



P(R = Y | V, hi) = Eh-! (P(R = Y| V, H)) 



X 7-1 / \ MN 

1 \ / 1 \ 




ex P (- £ f; fa-N-EL^M^^) 2 ^ (IU i 



\ w—l m—1 

which is again a multivariate Gaussian and can be expressed in 
closed form. Similarly, since P(R = Y | V, hi) is exponen- 
tial, J(V,W, hi) is also exponential and can be calculated 
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explicitly. However, both expressions are long and we do not 
give them here. 



